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Abstract 

We reanalyze deep inelastic scattering data of BCDMS Collaboration by in- 
cluding proper cuts of ranges with large systematic errors. We perform also the 
fits of high statistic deep inelastic scattering data of BCDMS, ShAC, NM and 
BFP Collaborations taking the data separately and in combined way and find 
good agreement between these analyses. We extract the values of the QCD cou- 



(N 



pling constant a s (M§) up to NLO level. The fits of the combined data for the 
| nonsinglet part of the structure function F 2 predict the coupling constant value 

^ ■ ol s {M 2 z ) = 0.1174 ± 0.0007 (stat) ±0.0019 (syst) ±0.0010 (normalization). The fits 

CN ■ of the combined data for both: the nonsinglet part of F 2 and the singlet one, lead to 

the values a,(Af|) = 0.1177±0.0007 (stat) ±0.0021 (syst) ±0.0009 (normalization). 
Both above values are in very good agreement with each other. 

(N [ 

2 ■ 1 Introduction 

o 

The deep inelastic scattering (DIS) leptons on hadrons is the basical process to study 
the values of the parton distribution functions (PDF) which are universal (after choosing 
of factorization and renormalization schemes) and can be used in other processes. The 
accuracy of the present data for deep inelastic structure functions (SF) reached the level 
at which the Q 2 -dependence of logarithmic QCD-motivated and power-like ones may be 
studied separately (for a review, see the recent papers M and references therein). 

In the present paper we review the results of our analysis at the next-to-leading 
^ ■ (NLO) order Q of perturbative QCD for the most known DIS SF F 2 (x,Q 2 ) taking into 
account SLAC, NMC, BCDMS and BFP experimental data We stress the power- 

like effects, so-called twist-4 (i.e. ~ 1/Q 2 ) contributions. To our purposes we represent 
the SF F 2 (x, Q 2 ) as the contribution of the leading twist part F2 QCD (x, Q 2 ) described by 
perturbative QCD and the nonperturbative part (i.e. twist-four terms ~ 1/Q 2 ): 

F 2 (x, Q 2 ) ee Ft U (x, Q 2 ) = Ff CD (x, Q 2 ) (l + (1) 

The SF F2 QCD (x, Q 2 ) obeys the (leading twist) perturbative QCD dynamics includ- 
ing the target mass corrections (and coincides with F 2 to2 (x, Q 2 ) when the target mass 
corrections are withdrawn). 

The Eq.(|ID allows us to separate pure kinematical power corrections, i.e. so-called tar- 
get mass corrections, so that the function h±{x) corresponds to "dynamical" contribution 

1 The evaluation of al(Q 2 ) corrections to anomalous dimensions of Wilson operators, that will be done 
in nearest future by Vermaseren and his coauthors (see discussions in ||), gives a possibility to apply 
many modern programs to perform fits of data at next-next-to-leading order (NNLO) of perturbative 
theory (see detail discussions in S). 
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of the twist-four operators. The parameterization ([I]) implies Qthat the anomalous dimen- 
sions of the twist-two and twist-four operators are equal to each other, that is not correct 
in principle. Meanwhile, in view of limited precision of the data, the approximation (Jl]) 
and one in the footnote 2 give rather good predictions (see discussions in [§]). 

Contrary to standard fits (see, for example, ||- [[K]]) when the direct numerical calcu- 
lations based on Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [llj are 
used to evaluate structure functions, we use the exact solution of DGLAP equation for 
the Mellin moments ilC 2 (Q 2 ) of SF F 2 tw2 (x,Q 2 ): 

M*{Q 2 )= [ x n ~ 2 F 2 k {x,Q 2 )dx (hereafter k = full,pQCD,tw2, ...) (2) 
Jo 

and the subsequent reproduction of Fl uU (x,Q 2 ), Ff (x,Q 2 ) and/or F|"' 2 (x, Q 2 ) at 
every needed Q 2 -value with help of the Jacobi Polynomial expansion method [13, O] (see 



also similar analyses at the NLO level |l4j] and at the NNLO level and above |15 



In this paper we consider in detail only fits of the nonsinglet part (NS) of SF F 2 and 
review only final results in general case, that can be found in 0. Moreover, we do not 
present exact formulae of Q 2 -dependence of SF F 2 which are also given in 0. We note 
only that the moments M^s{n,Q 2 ) at some Ql is theoretical input of our analysis and 
the twist-four term h^x) is considered as a set free parameters at each X{ bin. The set 
has the form h{ ree (x) = Yli=i h^Xi), where I is the number of bins. The constants hi{xi) 
(one per x-hm) parameterize x-dependence of h{ ree (x). 

2 Fits of F 2 : procedure 



Having the QCD expressions for the Mellin moments M^(Q 2 ) we can reconstruct the SF 
F 2 fc (x,Q 2 )as 

Afi max Tl 

F*> N ™(x,Q 2 ) = x a {l-x) b £ 9f (*) J>f (a,(3)Mj+2 {Q 2 ) , (3) 

71=0 j=0 

where 9^ ,b are the Jacobi polynomials f\ and a, b are the parameters, fitted by the condition 
of the requirement of the minimization of the error of the reconstruction of the structure 



functions (see Ref.|T3| for details). 

First of all, we choose the cut Q 2 > 1 GeV 2 in all our studies. For Q 2 < 1 GeV 2 , the 
applicability of twist expansion is very questionable. 

Secondly, we choose quite large values of the normalization point Qq. There are several 
reasons of this choice: 

• Our perturbative formulae should be applicable at the value of Qq. Moreover, the 
higher order corrections ~ «™(Qo) ( n — 2)> coming from normalization conditions 
of PDF, are less important at higher Qq values. 



2 The r.h.s. of the Eq.(|l|) is represented sometimes as F^ 00 (x,Q 2 ) + hi{x)/Q 2 . It implies that the 
anomalous dimensions of the twist-four operators are equal to zero. 

3 We would like to note here that there is similar method Jl6| to reproduce of structure functions, 
based on Bernstein polynomials. The method has been used in several analyses at the NLO level in JItJ 
and at the NNLO level in Ol. 
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It is necessary to cross heavy quark thresholds less number of time to reach Q 2 = 
Mf , the point of QCD coupling constant normalization. 

It is better to have the value of around the middle point of logarithmical range 
of considered Q 2 values. Then at the case the higher order corrections ~ (a s (Q 2 ) — 
Q-s{Ql)) n {n > 2) are less important. 



We use MINUIT program [ll| for minimization of two x 2 values: 



X 2 (F 2 



p ex P _ pteor i 2 



AF exp 



We would like to apply the following procedure: we study the dependence of x 2 / DOF 
value on value of Q 2 cuts for various sets of experimental data. The study will be done 
for the both cases: including higher twists corrections and without them. 

We use free normalizations of data for different experiments. For the reference, we use 
the most stable deuterium BCDMS data at the value of energy Eq = 200 GeV 4 . Using 
other types of data as reference gives negligible changes in our results. The usage of fixed 
normalization for all data leads to fits with a bit worser x 2 ■ 

3 Results of fits 

Hereafter we choose Q\ = 90 GeV 2 that is in good agreement with above conditions. We 
use also N max = 8, the cut 0.25 < x < 0.8, where the nonsinglet evolution is dominant. 

3.1 BCDMS C 12 + H 2 + D 2 data 

We start our analysis with the most precise experimental data obtained by BCDMS 
muon scattering experiment at the high Q 2 values. The full set of data is 607 points. 

It is well known that the original analyses given by BCDMS Collaboration itself (see 
also Ref. f9j) lead to quite small values a s (M|) = 0.113. Although in some recent papers 
(see, for example, || |2D|) more higher values of a s (M|) have been observed, we think 
that an additional reanalysis of BCDMS data should be very useful. 



Based on study (see also p0]| ) we proposed in [Qj that the reason for small values 
of a s (M§) coming from BCDMS data was the existence of the subset of the data having 
large systematic errors. We studied this subject by introducing several so-called F-cuts 
(see |^1|). Excluding this set of data with large systematic errors leads to essentially 



larger values of a s (M|) and very slow dependence of the values on the concrete choice of 
the F-cut (see below). 

We studied influence of the experimental systematic errors on the results of the QCD 
analysis as a function of Y cut3 , Y cuti and Y cut5 applied to the data. We use the following 
x-dependent y-cuts: 

y > 0.14 when 0.3 < x < 0.4 
y > 0.16 when 0.4 < x < 0.5 



4 Eq is the initial energy lepton beam. 

5 Hereafter we use the kinematical variable Y — (Eo — E)/Eo, where Eq and E are initial and scattering 
energies of lepton, respectively. 
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Figure 1: The study of systimatics at different Y cut values. The QCD analysis of BCDMS 
C 12 , H2, D2 data: the inner (outer) error-bars show statistical (systematic) errors. 




V > Ycut3 when 0.5 < x < 0.6 

V > YcuU when 0.6 < x < 0.7 

y>Y cut5 when 0.7<x<0.8 (4) 
and several sets iV of the values for the cuts at 0.5 < x < 0.8 given in the Table. 
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Table. The values of Y cut3 , Y cut4 and Y cut5 . 

The systematic errors for BCDMS data are given || as multiplicative factors to be 
applied to F 2 (x, Q 2 ): f r , fb, fs, fd and fh are the uncertainties due to spectrometer res- 
olution, beam momentum, calibration, spectrometer magnetic field calibration, detector 
inefficiencies and energy normalization, respectively. 

For this study each experimental point of the undistorted set was multiplied by a 
factor characterizing a given type of uncertainties and a new (distorted) data set was 
fitted again in agreement with our procedure considered in the previous section. The 
factors (f r , fb, f s , fa, fh) were taken from papers (see CERN preprint versions in 0). 
The absolute differences between the values of a s for the distorted and undistorted sets of 
data are given in the Fig. 1 as the total systematic error of a s estimated in quadratures. 
The number of the experimental points and the value of a s for the undistorted set of F2 
are also presented in the Fig. 1. 

From the Fig. 1 we can see that the a s values are obtained for N = 1 4- 6 of Y cut s, 
Ycuti and Y cut5 are very stable and statistically consistent. The case N = 6 reduces the 
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systematic error in a s by factor 1.8 and increases the value of a s , while increasing the 
statistical error on the 30%. 

After the cuts have been implemented (hereafter we use the set N = 6), we have 452 
points in the analysis. Fitting them in agreement with the same procedure considered in 
the previous Section, we obtain the following results: 

a s (M§) = 0.1153 ± 0.0013 (stat) ± 0.0022 (syst) ± 0.0012 (norm), 
= 0.1153 ± 0.0028 (total experimental error) 

where hereafter the symbol "norm" marks the error of normalization of experimental 
data. The total experimental error is squared root of sum of squares of statistical error, 
systematic one and error of normalization. 

3.2 SLAC, BCDMS, NMC and BFP data 

After these cuts have been incorporated (with N = 6) for BCDMS data, the full set of 
combine data is 797 points. 

To verify the range of applicability of perturbative QCD, we analyze firstly the data 
without a contribution of twist-four terms, i.e. when F 2 = . We do several fits 

using the cut Q 2 > Q 2 cut and increase the value Q 2 cut step by step. We observe good 
agreement of the fits with the data when Q 2 ut > 10 GeV 2 (see the Fig. 2). 

Later we add the twist-four corrections and fit the data with the usual cut Q 2 > 1 
GeV 2 . We have find very good agreement with the data. Moreover the predictions for 
a s (M z ) in both above procedures are very similar (see the Fig. 2). 

So, the analysis of combine SLAC, NMC, BCDMS and BFP data are given the fol- 
lowing results: 

• When twist-four corrections are not included and the cut of Q 2 is 10 GeV 2 at the 
free normalization 

X 2 /DOF = 0.98 and a s (M 2 z ) = 0.1170 ± 0.0009 (stat) (5) 

• When twist-four corrections are included and the cut of Q 2 is 1 GeV 2 

X 2 /DOF = 0.97 and a s {M 2 z ) = 0.1174 ± 0.0010 (stat) (6) 

Thus, as it follows from nonsinglet fits of experimental data, perturbative QCD works 
rather well at Q 2 > 10 GeV 2 . 

4 Summary 

We have demonstrated several steps of our study || of the Q 2 -evolution of DIS structure 
function F 2 fitting all modern experimental data at Bjorken variable x values: x > 10~ 2 . 

From the fits we have obtained the value of the normalization a s (M z ) of QCD coupling 
constant. First of all, we have reanalyzed the BCDMS data cutting the range with large 
systematic errors. As it is possible to see in the Fig. 1, the value of a s (M z ) rises strongly 
when the cuts of systematics were incorporated. In another side, the value of a s (M z ) 
does not dependent on the concrete type of the cut within modern statistical errors. 

We have found that at Q 2 > 10-^15 GeV 2 the formulae of pure perturbative QCD (i.e. 
twist-two approximation together with target mass corrections) are in good agreement 
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Figure 2: The values of a s (M|) and \ 2 & t different Q 2 -values of data cutes. The 
black points show the analyses of data without twist-four contributions. The white point 
corresponds to the case where twist-four contributions were added. Only statistical errors 
are shown. 

with all data. [] The results for a s (M z ) are very similar (see ||) for the both types of 
analyses: ones, based on nonsinglet evolution, and ones, based on combined singlet and 
nonsinglet evolution. They have the following form: 

• from fits, based on nonsinglet evolution: 

a s (M 2 z ) = 0.1170 ± 0.0009 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (7) 

• from fits, based on combined singlet and nonsinglet evolution: 

a s (M§) = 0.1180 ± 0.0013 (stat) ± 0.0021 (syst) ± 0.0009 (norm), (8) 

6 We note that at small x values, the perturbative QCD works well starting with Q 2 = 1.5-^2 GeV 2 and 
higher twist corrections are important only at very low Q 2 : Q 2 ~ 0.5 GeV 2 (see |2^| and references 
therein). As it is was observed in |24|, ^5) (see also discussions in [^2[ |23|, ^6|) the good agreement between 
perturbative QCD and experiment seems connect with large effective argument of coupling constant at 
low x range. 
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When we have added twist-four corrections, we have very good agreement between 
QCD (i.e. first two coefficients of Wilson expansion) and data starting already with 
Q 2 = 1 GeV 2 , where the Wilson expansion should begin to be applicable. The results for 
a s (M|) coincide for the both types of analyses: ones, based on nonsinglet evolution, and 
ones, based on combined singlet and nonsinglet evolution. They have the following form: 

• from fits, based on nonsinglet evolution: 

a a (M%) = 0.1174 ± 0.0007 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (9) 



• from fits, based on combined singlet and nonsinglet evolution: 

a s (M 2 z ) = 0.1177 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0009 (norm), (10) 

Thus, there is very good agreement (see Eqs. (0), @ and ([H])) between results 
based on pure perturbative QCD at quite large Q 2 values (i.e. at Q 2 > 10-^15 GeV 2 ) and 
the results based on first two twist terms of Wilson expansion (at Q 2 > 1 GeV 2 , where 
the Wilson expansion should be applicable). 

We would like to note that we have good agreement also with the analysis ]20| of 
combined HI and BCDMS data, which has been given by HI Collaboration very recently. 
Our results for a s (M|) are in good agreement also with the average value for coupling 



constant, presented in the recent studies (see || £7], [18], |28fl and references therein) and 
in famous Altarelli and Bethke reviews 
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